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Current theories of Kondo insulators employ the interaction of conduction electrons with localized 
Kramers doublets originating from a tetragonal crystalline environment, yet all Kondo insulators 
are cubic. Here we develop a theory of cubic topological Kondo insulators involving the interaction 
of Ts spin quartets with a conduction sea. The spin quartets greatly increase the potential for 
strong topological insulators, entirely eliminating the weak-topological phases from the diagram. 
We show that the relevant topological behavior in cubic Kondo insulators can only reside at the 
lower symmetry X or M points in the Brillouin zone, leading to three Dirac cones with heavy 
quasiparticles. 
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Our classical understanding of order in matter is built 
around Landau's concept of an order parameter. The 
past few years have seen a profound growth of interest in 
topological phases of matter, epitomized by the quantum 
Hall effect and topological band insulators, in which the 
underlying order derives from the non-trivial connected- 
ness of the quantum wave-function, often driven by the 
presence of strong spin-orbit coupling [IH5] . 

One of the interesting new entries to the world of 
topological insulators, is the class of heavy fcrmion, or 
"Kondo insulators" [TOHHj] . The strong-spin orbit cou- 
pling and highly renormalized narrow bands in these in- 
termetallic materials inspired the prediction [12] that a 
subset of the family of Kondo insulators will be Z2 topo- 
logical insulators. In particular, the oldest known Kondo 
insulator SmBg 17 with marked mixed valence character, 
was identified as a particularly promising candidate for a 
strong topological insulator (STI): a conclusion that has 
since also been supported by band-theory calculations 
p~3l [T6] . Recent experiments [T8H20] on SmBg have con- 
firmed the presence of robust conducting surfaces, large 
bulk resistivity and a chemical potential that clearly lies 
in the gap providing strong support for the initial pre- 
diction. 

However, despite these developments, there are still 
many aspects of the physics in these materials that are 
poorly understood. One of the simplifying assump- 
tions of the original theory [T^] was to treat the /- 
states as Kramer's doublets in a tetragonal environ- 
ment. In fact, the tetragonal theory predicts that strong 
topological insulating behavior requires large deviations 
from integral valence, while in practice Kondo insulators 
are much closer to integral valence [TT]. Moreover, all 
known Kondo insulators have cubic symmetry, and this 
higher symmetry appears to play a vital role, for all ap- 
parent "Kondo insulators" of lower symmetry, such as 
CeNiSn[2T] or CeRuASn fi [22] have proven, on improving 
sample quality, to be semi-metals. One of the important 



effects of high symmetry is the stablization of magnetic f- 
quartets. Moreover, Raman[23 j experiments and various 
band-theory studies [Ml [35] that it is the Kondo screening 
of the magnetic quartets that gives rise to the emergence 
of the insulating state. 
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FIG. 1: Contrasting the phase diagram of tetragonal [26] and 
cubic topological Kondo insulators. Cubic symmetry extends 
the STI phase into the Kondo limit. For SmB6 v = 3 — rif 
gives the valence of the Sm ion, while n / measures the number 
of /-holes in the filled 4/ 6 state, so that n/ = 1 corresponds 
to the 4/ 5 configuration. 

Motivated by this observation, here we formulate a 
theory of cubic topological Kondo insulators, based on 
a lattice of magnetic quartets. We show that the pres- 
ence of a spin-quartet greatly increases the possibility of 
strong topological insulators while eliminating the weak- 
topological insulators from the phase diagram, Fig. [T] 
We predict that the relevant topological behavior in sim- 
ple cubic Kondo insulators can only reside at the lower 
point group symmetry X and M points in the Brillouin 
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zone (BZ), leading to a three heavy Dirac cones at the 
surface. One of the additional consequences of the un- 
derlying Kondo physics, is that the coherence length of 
the surface states is expected to be very small, of order 
a lattice spacing. 

While we outline our model of cubic Kondo insulators 
with a particular focus on SmB 6 , the methodology gen- 
eralizes to other cubic Kondo insulators. SmB 6 has a 
simple cubic structure, with the B 6 clusters located at 
the center of the unit cell, acting as spacers which medi- 
ate electron hopping between Sm sites. Band-theory [25] 
and XPS studies [27] show that the 4f orbitals hybridize 
with d-bands which form electron pockets around the X 
points. In a cubic environment, the J = 5/2 orbitals split 
into a Tj doublet and a Tg quartet, while the fivefold de- 
generate d-orbitals are split into double degenerate e g 
and triply degenerate <2g orbitals. Band theory and Ra- 
man spectroscopy studies [53] indicate that the physics 
of the 4/ orbitals is governed by valence fluctuations in- 
volving electrons of the Tg quartet and the conduction 



e g states, e" + 4/ 5 (r 8 Q) ) ^ 4/ 6 . The Tf {a = 1,2) 
quartet consists of the following combination of orbitals: 
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|4 2) ) = |±|>. This then 
leads to a simple physical picture in which the Tg quartet 
of /-states hybridizes with an e g quartet (Kramers plus 
orbital degeneracy) of d-states to form a Kondo insulator. 

To gain insight into how the cubic topological Kondo 
insulator emerges it is instructive to consider a simplified 
one-dimensional model consisting of a quartet of con- 
duction d-bands hybridized with a quartet of /-bands 
(Fig. [2^,) . In one dimension there are two high symmetry 
points: F (k = 0) and X (k — tt), where the hybridiza- 
tion vanishes [12 HH [5S]). Away from the zone center 
r, the /— and d— quartets split into Kramers doublets. 
The Z2 topological invariant v\r> is then determined by 
the product of the parities V\t> — Sr^x of the occupied 
states at the T and X points. However, the /-quartet 
at the r point is equivalent to two Kramers doublets, 
which means that <5r = (±1) 2 is always positive, so that 
v\r) = 5x and a one-dimensional topological insulator 
only develops when the / and d bands invert at the X 
point. 

Generalizing this argument to three dimensions we see 
that there are now four high symmetry points T, X, M 
and R. The /-bands are fourfold degenerate at both T 
and R points which guarantees that 5r — Sr = +1 
(Fig. [2}}). Therefore, we see that the 3D topological 
invariant is determined by band inversions at X or M 
points only, v?,d = (SxSm) 3 = SxSm ■ If there is a band 
inversion at the X point, we get z/ 3 £> = SxSm = — 1- 
In this way the cubic character of the Kondo insulator 
and, specifically, the fourfold degeneracy of the /-orbital 
multiplet protects the formation of a strong topological 
insulator. 

We now formulate our model for cubic topological 




FIG. 2: Schematic band structure illustrating (a) ID Kondo 
insulator with local cubic symmetry and (b) 3D cubic Kondo 
insulator. Hybridization between a quartet of d-bands with a 
quartet of f-bands leads to a Kondo insulator. The fourfold 
degeneracy of the /- and d-bands at the high symmetry T 
and R points of the Brillouin zone guaranties that the 3D 
topological invariant is determined by the band inversions at 
the X and M points only. 



Kondo insulators. At each site, the quartet of f and d 
-holes is described by an orbital and spin index, denoted 
by the combination A = (a, a) (a = 1,2, a = ±1). The 
fields are then given by the eight component spinor 



Xox(j) 



(1) 



where d\(j) destroys an d-hole at site j, while X \(j) = 
|4/ 6 )(4/ 5 ,A| is the the Hubbard operator that destroys 
an /-hole at site j. The tight-binding Hamiltonian de- 
scribing the hybridized f-d system is then 



H = X>l(i)/iAv(Ri - R,-)*v0") 



(2) 



in which the nearest hopping matrix has the structure 



h(R) 



fh d (R) V(R) 
l0(R) hf(R) 



(3) 



where the diagonal elements describe hopping within the 
d- and /- quartets while the off-diagonal parts describe 
the hybridization between them, while R £ (±x, ±y, ±z) 
is the vector linking nearest neighbors. The various ma- 
trix elements simplify for hopping along the z-axis, where 
they become orbitally and spin diagonal: 



h l (z) = t l 1 ' 



m 



V(z) = iV 







(4) 



where I = d, f and rji is the ratio of orbital hopping 
elements. In the above, the overlap between the Tg 1 -* 
orbitals, which extend perpendicular to the z-axis is ne- 
glected, since the hybridization is dominated by the over- 
lap of the the Tg ; orbitals, which extend out along the 
z-axis. The hopping matrix elements in the x and y 
directions are then obtained by rotations in orbital/spin 
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space, so that h(x) — U y h(z)U^ and h(y) — U- x h(z)U^_ x 
where U y and U- x denote 90° rotations about the y and 
negative x axes, respectively. 

The Fourier transformed hopping matrices h(\s) = 
J2r h(R)e~ lk ' IL can then be written in the compact form 

where I — d, f. Here e l are the bare energies of the 
isolated d and f-quartets, while <fii (k) = 



(k) = 3(c x + c y ) and </> 3 (k) = V3(c 



-y 



(c a cos k a ,a — x,y, z). The hybridization is given by 



V(k) = 



If 3(<J X +i<Jy) ^/i((T x - lljy) 

6 \V3(a x - id y ) a x + ia y + 4ct 2 



(6) 



where we denote a a — cr a sin k a . Note how the hy- 
bridization between the even parity d-states and odd- 
parity f-states is an odd parity function of momentum 

v(k) = -y(-k). 

To analyze the properties of the Kondo insulator, we 
use a slave boson formulation of the Hubbard operators, 
writing X\ (j) = f{(j)b 3 , where /jjO) = |4/ 5 ,A) creates 
an f-hole in the T 8 quartet while &t j 0) = |4/ 6 ) denotes 
the singlet filled 4/ shell, subject to the constraint Qj — 
b ] b J+Y,\ fjxfjX = 1 at each site. 
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FIG. 3: Temperature dependence of the hybridization gap 
parameter b and the renormalized /-level position (inset) for 
various values of the bare hybridization (see Supplementary 
Materials for more details). 

We now analyze the properties of the cubic Kondo in- 
sulator, using a mean-field treatment of the slave boson 
field bi, replacing the slave-boson operator hi at each site 
by its expectation value: (bi) — b so that the f- hopping 
and hybridization amplitude are renormalized: tf — > b 2 tf 
and Vdf —> bVdf- The mean- field theory is carried out, 
enforcing the constraint b 2 + (nj) = 1 on the average. 
In addition, the chemical potentials td and £/ for both 
d-electrons and /-holes are adjusted self-consistently to 



produce a band insulator, rid + nj = 4. This condi- 
tion guarantees that four out of eight doubly degenerate 
bands will be fully occupied. The details of our mean- 
field calculation are given in the Supplementary Materials 
section. Here we provide the final results of our calcula- 
tions. 

In Fig. [3] we show that the magnitude b reduces with 
temperature, corresponding to a gradual rise in the Sm 
valence, due to the weaker renormalization of the /- 
electron level. The degree of mixed valence of Sm + is 
given then by v = 3 — (rif). In our simplified mean- 
field calculation, the smooth temperature cross-over from 
Kondo insulating behavior to local moment metal at 
high temperatures is crudely approximated by an abrupt 
second-order phase transition. 




FIG. 4: Band structure consistent with PES and LDA studies 
of SmBe computed with the following parameters: nj = 0.48 
(or b = 0.73), V = 0.05 eV,t d = 2 eV, /i d = 0.2 eV , rj = 
rf = -.3, e f = -0.01 eV (e /0 = -0.17 eV), t f = -.05 eV, 
T = 10~ 4 eV and the gap is A = 12 meV. Shaded region 
denotes filled bands. Inset shows the ground-state energy 
computed for a slab of 80 layers to illustrate the three gapless 
surface Dirac excitations at the symmetry points F, X'andX" . 



Fig. |4]shows the computed band structure for the cubic 
Kondo insulator obtained from mean-field theory, show- 
ing the band inversion between the d- and /-bands at 
the X points that generates the strong topological insu- 
lator. Moreover, as the value of the bare hybridization 
increases, there is a maximum value beyond which the 
bands no longer invert and the Kondo insulator becomes 
a conventional band insulator. 

One of the interesting questions raised by this work 
concerns the many body character of the Dirac electrons 
on the surface. Like the low-lying excitations in the va- 
lence and conduction band, the surface states of a TKI in- 
volve heavy quasiparticles of predominantly /-character. 
The characteristic Fermi velocity of these excitations 
v* F = Zvf is renormalized with respect to the conduc- 
tion electron band group velocities, where Z — m/m* is 
the mass renormalization of the f-electrons. In a band 
topological insulator, the penetration depth of the sur- 
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face states £ ~ Vf/A, where A is the band-gap, scale 
that is significantly larger than a unit-cell size. Paradox- 
ically, even though the Fermi velocity of the Dirac cones 
in a TKI is very low, we expect the characteristic pene- 
tration depth £ of the heavy wavefunctions into the bulk 
to be of order the lattice spacing a. To see this, we note 
that £ ~ -yfi- , where the indirect gap of the Kondo insula- 
tor A g is of order the Kondo temperature A fl ~ Tk . But 
since Tk ~ ZW , where W is the width of the conduction 
electron band, this implies that the penetration depth of 
the surface excitations £ ~ vf/W ~ a is given by the 
size a of the unit cell. Physically, we can interpret the 
surface Dirac cones as a result of broken Kondo singlets, 
whose spatial extent is of order a lattice spacing. This 
feature is likely to make the surface states rather robust 
against the purity of the bulk. 

Various interesting questions are raised by our study. 
Conventional Kondo insulators are most naturally un- 
derstood as a strong-coupling limit of the Kondo lat- 
tice, where local singlets form between a commensurate 
number of conduction electrons and localized moments. 
What then is the appropriate strong coupling description 
of topological Kondo insulators, and can we understand 
the surface states in terms of broken Kondo singlets? A 
second question concerns the temperature dependence of 
the hybridization gap. Experimentally, the hybridization 
gap observed in Raman studies [23J is seen to develop in 
a fashion strongly reminiscent of the mean-field theory. 
Could this indicate that fluctuations about mean-field 
theory are weaker in a fully gapped Kondo lattice than 
in its metallic counterpart? 

We end with a few comments on the experimental con- 
sequences of the above picture. One of the most dramatic 
consequences of the quartet model is the prediction of 
three Dirac cones of surface excitation of substantially 
enhanced effective mass. By contrast, were the underly- 
ing ground-state of the f-state a Kramers doublet, then 
we would expect a single Dirac cone excitation. These 
surface modes should be observable via various low en- 
ergy spectroscopies. For instance, the heavy mass of the 
quasi particles should appear in Shubnikov-de Haas or cy- 
clotron resonance measurements. The quasiparticle mass 
to* = hkp/vp ~ m e {kFa)Z~ l will depend on the Fermi 
momentum kp of the Dirac cones. Scanning tunneling 
spectroscopy measurements of the quasiparticle interfer- 
ence created by the Dirac conese, and high resolution 
ARPES measurements may provide a direct way to ob- 
serve the predicted three Dirac cones. 

To summarize, we have studied the cubic topological 
Kondo insulator, incorporating the effect of a fourfold de- 
generate /-multiplet. There are two main effects of the 
quartet states: first, they allow the low fractional filling 
of the band required for strong topological insulating be- 
havior to occur in the almost integral valent environment 
of the Kondo insulator; second, they double the degen- 



eracy of the band-states at the high-symmetry T and R 
points in the Brillouin zone, removing these points from 
the calculation of the Z2 topological invariant so that the 
only important band crossings occur at the the X or M 
points. This leads to a prediction that three heavy Dirac 
cones will form on the surfaces [T3J [TB] . 
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tute Award, Grant DMR-0844115. 
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SUPPLEMENTARY MATERIALS FOR CUBIC TOPOLOGICAL KONDO INSULATORS 

These notes provide: 

• details of the derivation of the tight-binding Hamiltonian for a cubic Kondo insulator. 

• derivation of the mean-field theory for the infinite U limit 

• derivation of the mean-field equations. 

Rotation matrices 

To construct the Hamiltonian, we evaluate the hopping matrices along the z-axis, and then carry out a unitary 
transformation to evaluate the corresponding quantities for hopping along the x and y axes. 
Consider a general rotation operator 

K = RA, (7) 

where R and A describe tt/2 rotations about a principle axis of the crystal in real and spin space respectively. The 
Hamiltonian in a cubic environment is invariant under these transformations: H = TZ^HTZ. We now write the 
directional dependence of the Hamiltonian explicitly. 

H(x, y, z) = (RA)^H(x, y, z)RA = A\R^H(x, y, z)R)A. (8) 

We can always choose the rotation 1Z to transform in the cyclic manner: x — > y — > z — > x, then substituting z = y = 0, 

H(x, 0, 0) = A\RlH(x, 0, 0)R x )A = A* #(0, 0, x)A, (9) 

hence we can assume the 3D Hamiltonian of the following form 

H 3D = H(k z ) + AH(k x )tf + A^H(k y )A, (10) 

Where to rotate in the opposite direction we use A -1 = A*. Using the Wigner D-functions we can construct the 
rotation in angular momentum space 

A = e -i<xJ* e -iPJy e -ilJ*, (11) 

Here, A denotes the rotation operator in the "ZYZ" convention, corresponding to a rotation around the z-axis, followed 
by a rotation around the y-axis, and then the new - z-axis. The matrix elements of this operator are then 

D J m , m = (Jm'\A\Jm) 

= [(j + - m)!(j + m')!(j - m')!] 1/2 

x v tiZ 

{j-m 1 -x)!0'+m-x)!(x+m' -m)!x! 

2j+m—m.' — 2x / a \m'—m+2x 

X I COS " 



(cosf)'"* (-sin I)' 



so that upon rotating the state \Jm) using R x we obtain A^= 5 /2) = D^,(0,7r/2,7r/2), and A^= 2 ) = D 2 mm , ® 

1 /2 

D ' i (0, 7r/2,7r/2). One can now obtain the transformation matrices for a given multiplct (e g ,i 2 g,r 8 etc). For the e g 
doublet we are to read off the matrices M^- = (i\A x \j) from the transformation equation 

A x \e g :m) = \e g :n)M* m , (12) 

where \e g : m) = {\d x 2_ y 2 f), \d z 2 f), |d a; 2_ 1/ 2 J,), |rf z 2 |)} is an e ff doublet. Similarly for the r 8 quartet, 

A x \T s ,a) = \T 8 ,p)M f 0a , (13) 
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where the quartet is denoted by |rg,a), where a £ [1,4]. The result is, 

/ _J_ Ij.fs L_ _i,* 4 /i"\ / _J_ 



M a mn = e 



i I,' /3 t_ _1- /?\ 

2% /2 2'\/2 2V2 2'Y2 

2 l \/ 2 2V2 2'Y 2 2^2 



V 



1 V 2 

2v^2_ 2 
2 2V2 



1 

2\/2 



1 

2\/2 



2V2 


2 


2V2 2 




1 


VI 1 


2 


2V2 


2 2\/2 


i 






2V2 


2 l \l 2 


2\/2 2 4 \/2 




i 


2*V 2 2^ / 


2 i y 2 


2V2 



(14) 



Note that the rotation cyclically exchanges x — >■ y z x, thus applying it three times gives overall " — 1" due to 
fermionic statistics. 

M 3 = -1. 



Finally, we can use (14) to derive the hamiltonian defined in (5) and (6) of the main paper. 

Details of the mean-field theory 

The full Hamiltonian for the problem also contains a term describing the local Hubbard interaction between the 
/-electrons: 



i a — 1 /3^a 



(15) 



We consider the infinite U , where we can project out all states with occupation number larger than one by replacing the 
bare f-electron fields by Hubbard operators. We represent the Hubbard operators using a slave boson representation, 
as follows: 



fL -> X a0 (i) = fjji, f ia -> Xoa(i) = b\f ia , 

supplemented by a constraint of no more than one /-electron per each site (U = 00): 



(16) 



E flf*° + b ^ = L 

a=l 

The partition function corresponding to the model Hamiltonian H = H c + H t + Hhyb above and with constraint 

7T//3 



condition (17) reads : 



f3d\ 
n 



P(6,6t ) /,/t )C)C t) exp _ / L {r)dT 



(18) 



where the Lagrangian L(r) is 

4 



L =E^+E E +m^&}U /j +x; E 4.(Jr+^)W 

i ij a,/3=l L V / J kcr a b=1 V / 

1 2 4 — > / 4 \ 

+ 2^Y,J2J2 (Vtojpcirbtfjf, + h.c.) + e e 4/ia + b h 1 

(ij) kcr 0=1/9=1 j \a=l / 

Mean-field (saddle-point) approximation corresponds to the following values of the bosonic fields: 

kq( r ) = W q,o, *A q (r) = (Ef - £/)6 q ,o, 



(19) 



(20) 



where both a and ej are r- independent. In the mean-field theory we choose a "radial" gauge where the the phase 
of the &-field has been absorbed into the f-electron fields. Also, for an insulator, we need a filled quartet of states at 
each site, so that 



n c + rif = 4. 



(21) 



Note, the parameter b also renormalizes the /-hopping elements. Indeed, it follows that when i ^ j, fj a fjp — > 
X a o(i)Xof;(j) — fi a bitfjfj[). However, the onsite occupancy is unrenormalized by the slave boson fields, since in the 

infinite U limit, the onsite occupancy X aa (i) = fj a fi a - 

The first three terms together with the last term in the Lagrangian ( 19 ) can be written using the new fermionic 
basis and (20 1. It follows: 



Lo{r) = (6 2 - l)(E f - e f ) + /W * + E%> f nka 



kcr a=l 



dr 



E 



id) , , 
k I u akcr 



where the renormalized /-electron dispersion is 



E 



(/) 



E f + 2t f b 2 (c x + c y + c z + (-1)™^/ c| + c% + c| - c x Cy - CyC z - c z c x \ , n = 1,2. 



(22) 



(23) 



Correspondingly, the hybridization matrix in the new fermionic basis can be obtained via unitary transformation with 
the matrix 



I M k Wk \ 

itk «k 

— «k Mk 

\ -u k u k / 



/uk -«k \ 

Uk — Uk 

«k Uk 

\0 Dk u k J 



(24) 



Thus we need to calculate the elements of the matrix 



Hv = iVdf-Y, d\U^H hvb 0<.)U k f k + h.c. 



(25) 



We can use the following relations 



«k - «k = 



2-Rk 



V3 

"k«k = 7^r( c z ~ c v)- 



(26) 



In our subsequent discussion it will be convenient to write Hy (k) in a more compact form. To derive the corresponding 
expression we first recall that Hh y b can be written as follows 



H hyb (k) =iVdf~ ^ x {k)®a x + 4> v (k)®a y + j> z (k)®a z , 4> a = <#*t + <pfr z + <^f x , 

-2 1 1 

4>q = -(sin k x , — sin k y , sin k z ), <f>\ = -(sin fc^, — sin k y , —2 sin fc z ), 02 = H — ^=(sin fc x , sin k y , 0). 



(27) 



Again, note the change from minus sign to plus sign in front of <p2 term. This yields the agreement with the 
hybridization Hamiltonian obtained from the rotations method. We also have 



£4: = Ukfo ® (To + iVkfy ® O"0, ^k 1 = "k^o <8> CTq — i«kT y <8> <7 , 



After some algebra we find 



W k 1 (f ® o-z) W k = r ® a x 



U k 1 (f 2 ® a x )U k = (ul - vl)(f z ® + 2u k Wk(T a; ® o^), 
^k 1 ( T * ® ct x )Z4 = (ul - vl)(r x ® <7 X ) - 2u k «k(T z «) <7 X ). 

These results become much more transparent if we express u k and v k in terms of the angle 9 k : 

u k = cos(6» k /2), w k = sin(6» k /2). 

Then we see that 



1 (>(k)®^)w k = <H 



k) <g> o- x , 



(28) 



(29) 



(30) 



(31) 
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where now 



cos (?k — sin (9k 
sin 0k cos 6*k 



Similarly, we find 



Wk 1 (>(k)®«7 a 
$ y = $»t + f z + 

Thus, we can re- write (251 as follows 



W k = $ a (k) ® fr„, U k 1 (>(k) ® a 2 ) W k = * z (k) ® a z 



s>y — eh y 



n = 4> z , 



cos ^k sin 
— sin 0k cos 6*k 

cos 6*k — sin 0k 
sin 0k cos #k 



H v = iV d A 5I4a[*k]a^/k/9 + h.c, $ k = **(k) ® <T X + $ s (k) ® & y + <l z (k) <g> a z 

k 

Next we derive the mean-field equations. 



(32) 



(33) 



(34) 



Derivation of the mean-field equations 

To derive the mean field equations we first integrate out (i-electrons by making the following change of variables in 
the path integral: 



4 4 + fl^Gd^, k), d k -> 4 + G d (zw, k)l»k/k 

(iu-E^ \ 



Q^k) 



V 



iuj — E 





(<0 

lk 







iuj — E. 



(d) 



2k 

iu> 



(35) 



Then the resulting action is Gaussian and the /-electrons can be integrated out. This yields the effective action of 
the form 



S eff = {b 2 - \){E f - £/ )-T^^logdet \Gj}(iu,,k)} , Gjj = df^k) - (V^-f^G^, k)$ k . (36) 



iuj k 



The renormalized /-electron correlation function Gjj has a block diagonal form: 



/ a-} 



G ff 



(V d fb) 2 
4 







A* 



G X f{iu,k) 
A 



A* 



(V d /b) 2 

A 2k 

u—hi^ 



V iw-.E& J i^ig 5 ^-B[ d k » 



iuj—E. 
A* 



777T 



ioj—E~ 



a; w 



iu — E 



(d) 



(Vd/fc) 2 



A* 



a;,. 



iu~E 



(d) 



iuj—E. 



(d) 



(^fc) 2 y 



(37) 



where the diagonal elements are given by 



er / i (z W ,k) = JW -^( ) -(y d/ 



g 2 7 1 (z W ,k) = lw -^{ ) -(v/ d/ ^) 





$2k 




iu - E y ^ 


*2k 


, v ik 


iu; - E lk 


iUJ ~ ^2k 



(38) 
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and we have introduced the following functions 



(±) 

lk 



(*g ± + m ± <s>xy + m ± $?) 2 , $ 2k - ($^) 2 + ($d 2 + ($^) 2 , 



<i> iiv 

A 2k - $f - $f) + ($| - - $f + - *?)], 

A lk - ($f - i$i)($g + $f) - $f[<&§ + $f - + $?)], 
A 2k - *f [Sg - $f - - - (*f - i$i)($g - $?). 



(39) 



Then the determinant of the matrix G// is 



dot 



gj* (iw, k) = \ g-j k)^-/ (»«, k) - - ) 



-(v df -y 



A 



lk 



2 y 
A 2k 



Au 

iui — E 



A 



(d) 

lk 

A* 

^lk 



2k 1 

A; k 



Aj t 

iui — E 



(d) 

lk 



Ajk 



(d) 

2k , 



(40) 



iui — E 



(d) 

lk 



Furthermore, we find that 



iui — E. 



|A lk | 2 + |A lk | 2 - $W$ 2k , 



(d) 

2k , 



iui — E 



|A 2k | 2 



i'-i) 

lk 



iui — E. 



id) 

2k , 



A 2k | 2 = $L } $ 



lk ^2k 



(41) 



so that certain terms in the expression (40) will cancel. Therefore, to derive the mean field equation we will need to 
find the saddle point of the following effective action: 

4 

S eff = (b 2 - l)(E f - Ef )-2TJ2Y1 MAk^; E,,b)], P 4 kM = ]J(iui - e 4k ). (42) 

iuj k 2—1 

where we have introduced the polynomial: 



P 4k M = (iui - E[Q){iu - E&>)(iui - E$>)(iu, - E%) + {V^-frt- 

- {V df h -f {[$ 2k (i w - E$) + (iw - E%)](iui Eg) + [<f 2k (zc - 
7k = + $ 2k - Ai k A; k - A* lk A 2k - A lk A* k - A^ k A 2k . 

Thus, two of the three mean field equations are formally given by 



^( ) ) + *ik(^-^ ) )](^-^k ) )}, 



(43) 



dS, 



eff 



dE f 



= 0, 



dS, 



eff 



db 



= 0. 



We have 



»=1 k 



EI ( £ ik ^ £;k) 



(44) 



(45) 



where the function /(e) is: 



/(A) = 2X 3 - \ 2 [2(E {d) 



lk 



2k J 



K2E®Ej* + (< + E%)(e£ + - (^) 2 ($ + + $- k + 2$ 2k ) 



Ad) 



(46) 



fc) 2 [($2k + + ($2k + - E$E$(El£ + 



■^ru^k + ^kj^k 

The derivation of the last mean-field equation can be compactly written as follows 

de 



8(E f ~e f ) + 2^2T,di n ^ = °- 



t=l k 



/ 



(47) 



